Laboratory experiments on the generation of internal tidal beams over steep slopes 
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We designed a simple laboratory experiment to study internal tides generation. We consider a 
steep continental shelf, for which the internal tide is shown to be emitted from the critical point, 
which is clearly amphidromic. We also discuss the dependence of the width of the emitted beam 
on the local curvature of topography and on viscosity. Finally we derive the form of the resulting 
internal tidal beam by drawing an analogy with an oscillating cylinder in a static fluid. 
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Over the past two decades it has become apparent that 
substantial internal tides can be generated by tidal cur- 
rents over ridges and other rough topography of the ocean 
floor. This problem is of paramount importance since 
baroclinic tides, generated by barotropic currents over 
ocean ridges and seamounts, are an important source of 
the ocean interior mixing. Indeed, recent observations 
suggest that mixing in the abyssal oceans is rather weak, 
except in localized regions near rough topography. It 
explains the discrepancy between the observed intensity 
of mixing in the interior of the oceans and is required 
to satisfy models of ocean circulation [7]. This suggests 
that these topographic effects must be incorporated into 
realistic climate and circulation models. 

This question has been studied theoretically and nu- 
merically, while related oceanic observations have been 
reported. Baines showed [2] first that points where the 
topographic slope coincides with the angle of propagation 
might be an effective generator of internal wave. Such a 
point is referred to as the critical point (see Fig. [T]). In- 
ternal beams emanating from the continental slope have 
been observed recently [J, H, [3, EE E3- 0n tne the- 
oretical side, let us mention several interesting descrip- 
tions [H, Ql| of abrupt discontinuity cases such as a to- 
pographic step, a knife-edge ridge (zero width ridge as a 
simple model for tall ridges such as the Hawaiian one), 
a tent-shaped ridges, or smooth ridges (gaussian, poly- 
nomial). Finally, recent numerical reports have al so p ro- 
vided some new insight about this mechanism Il0|. 

All these works converge to conclude that beams of 
internal tides energy arise due to the interplay between 
oscillating currents and bathymetric features. Of partic- 
ular interest is the location where the slope is critical, i.e. 
where the direction of the wave beam is tangent to the 
slope. In the remaining of the paper, we experimentally 
study the location of emission and the form of the in- 
ternal beam which is emitted. We also explain how this 
form can be predicted by drawing an analogy between 
the emission of an internal tide from a static topography, 
and the emission of internal waves by an oscillating cylin- 
der in a static fluid. The agreement is very good, even 
quantitatively. 



The experiment was performed in a rectangular- 
shaped tank (length 120 cm, height 40 cm and width 
10 cm) as sketched in Fig.[TJ A curved thin (0.5 cm) PVC 
plate, introduced before filling, played the role of the con- 
tinental shelf. The experiment was performed with a lin- 
ear stratification, pure water at the free surface while 
highly salted water at the bottom, resulting in a con- 
stant Brunt-Vaisala buoyancy frequency TV ~ 0.81 rad/s. 
In the absence of rotation, the dispersion relation of in- 
ternal gravity waves reads uo = N sin 0, where is the 
angle of the beam with respect to the horizontal. 




FIG. 1: Sketch of the experimental set-up. A PVC plate is os- 
cillating between the solid and dashed oblique lines indicated 
on the right, creating the tide. The continental shelf consists 
of two planar PVC plates connected by a quarter of cylinder 
of radius R = 3.3 cm. The upper plate is nearly horizontal, 
while the downgoing one makes an angle a = 78° with the 
horizontal. The aspect ratio h/H is 2/3, for a total water 
height h + H = 19 cm. The critical point C indicates the 
location where the internal tide generation is supposed to oc- 
cur, i.e. where the internal beam inclined at an angle — 23° 
with the horizontal is tangent to the topography. The dot- 
ted rectangle indicates the zone captured by the camera and 
presented in Fig. [2Jl. 

The tidal current was generated by a second PVC plate 
localized far from the continental slope, i.e. in the region 
corresponding to the abyssal plane or to the open ocean. 
This plate was fixed at its bottom end, with its top end 
oscillating horizontally. The oscillation Asm(ut) was 
monitored by a position control motor of a printer ap- 



paratus driven by a sinusoidal electronic signal delivered 
by a Kepco amplifier. The phase origin at t = corre- 
sponds to a vertical position of the paddle and the move- 
ment is initiated leftwards towards the topography. Sev- 
eral forcing tidal frequencies and amplitudes were tried 
out. The tidal frequency uo was satisfying the condition 
< uo < N, so that internal tides may be freely radiat- 
ing. Results and figures discussed throughout this paper 
correspond to uo/ (2ir) = 0.050 Hz for which = 23°, and 
A = 0.5 cm. 

The density gradient field within the tank was finally 
obtained with the standard synthetic Schlieren tech- 
nique [H by acquiring successive side views with a CMOS 
AVT Marlin F131B Camera. A sheet of paper with 
points randomly scattered was located 150 cm behind 
the tank. Successive images (1280 x 300 pixels) of those 
points obtained with the camera located 280 cm in front 
of the tank were adequately treated by a homemade soft- 
ware. A correlation image velocimetry algorithm [6[ was 
applied between a reference image taken prior to the ex- 
periment and the different snapshots. As we used 21x21 
pixels correlation boxes with a 75% overlap, the reso- 
lution for the two dimensional density gradient field is 
256 x 60. 

Figured show a snapshot of the vertical density gra- 
dient. To facilitate the visualization, a narrower domain 
than in the experiments is shown. It corresponds to the 
dotted rectangle depicted in Fig. [TJ To increase the visi- 
bility, we filtered at the excitation frequency uo the time 
series of the density gradient field over one experimen- 
tal tidal period (see Ref. [l2[ for more details about this 
method). 

An internal tide is clearly seen to emanate from the 
upper part of the slope and radiate away from the shelf- 
break in both directions. However, it is important to 
stress the absence of a third beam radiated transversally 
to the topography. This contradicts Baines analytical 
theory Q in which this third beam reflects on the surface 
and is present in the general solution for the downward 
propagating wave. The explanation presumably lie in the 
presence of a singular point in Baines' case or in Saint 
Laurent et al numerical simulation [l8j]. In numerical 
experiments with a smooth slope and a typically shallow 
continental shelf [9[, no such beam was found and the 
presence of this third beam was already ambiguous in 
previous experiments [3]. 

The amplitude of the vertical density gradient is 
given in terms of variations of the squared Brunt- 
Vaisala frequency. If one considers the original value of 
N 2 = 0.66 rad 2 /s 2 , the measured amplitude of AN 2 = 
±0.005 rad 2 /s 2 for the internal wave is a very small per- 
turbation of the original stratification. Using the mass 
conservation relation 

iuo(p - p) = w^- , (1) 

where p is the perturbated density, p(z) the initial density 
and w the vertical density, we obtain typical vertical ve- 
locities of order ±0.11 mm/s, corresponding to vertical 
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FIG. 2: (Color online) Panel a presents the two-dimensional 
vertical density gradient, while panel b shows the phase. Both 
plots have been averaged over one tidal period of the experi- 
ment around t — 3T. Horizontal and vertical distances are in 
cm. The generation point C appears to be an amphidromic 
point around which the phase of the wave rotates. 

displacements of ±0.15 mm. This is one of the inter- 
ests of the synthetic Schlieren technique that allows to 
measure very weak perturbations of the buoyancy field 
and thus to investigate weakly nonlinear regimes. This 
vertical amplitude has to be compared to the barotropic 
elevation of the water induced by the forcing. The paddle 
oscillates with an A = 0.5 cm amplitude, and the width 
of the free water volume in the tank is 90 cm for a height 
of H = 12 cm. The corresponding elevation of the water 
at the level of the slope is thus AH = 0.7 mm. The baro- 
clinic component observed is thus still a perturbation of 
the barotropic tide. 

To the right of the generation point C, energy propaga- 
tion is downward while, to the left, it is upward. As a con- 
sequence of the internal waves propagation law for which 
group and phase velocities are orthogonal with opposite 
vertical components, phase propagation is thus upward 
to the right and downward to the left. One can conclude 
from this simple observation that the phase has to rotate 
around C, which is therefore an amphidromic point. Our 
filtering technique allows to evaluate the phase of the 
wave [12], which is plotted in Fig. [2)3. One can clearly 
see that the isophase lines converge on a single point pre- 
viously referred as the generation point C, around which 
the phase rotates uniformally. 

Whereas the location and the inclination of the inter- 
nal tide are well understood, the selection mechanism of 
the width of the beam was not yet clearly identified by 
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previous studies. Several length scales can be consid- 
ered in this problem. The first one corresponds to the 
thickness of the oscillating boundary layer S = (yjuj) 1 ! 2 
where v is the kinematic viscosity. In the present case, 
5 = 1.8 mm. The second one is the local radius of cur- 
vature of the continental shelf R = 3.3 cm. At last, the 
dimensions of the shelf itself (ft, H...) that play a role Q 
in the "flat-bump" geometry for which a < happen to 
be irrelevant in the configuration of a steep topography. 

Our understanding is directly inspired from the gen- 
eration of internal waves by oscillating cylinders. One 
can therefore try to draw an analogy between the inter- 
nal tide generation by a curved static topography of a 
given radius of curvature R and the internal waves gen- 
eration by an oscillating cylinder of the same radius R in 
a stratified fluid. 



locity component v s for each beam is shown to be given 

by ' 




FIG. 3: (Color online) Comparison between experiment (a) 
and theory (b) for the vertical density gradient field. Panel (a) 
has been averaged over one tidal period around 3T. We use 
the Hurley and Keady analytical solution [l3[ for oscillating 
cylinders to modelize the internal wave beam. We consider 
the tangential circle at the shelf break as the generating cylin- 
der section and compute the analytical solution isolating the 
upper right emission point. Vertical dashed lines correspond 
to the cross section plotted in Fig. |4] 

In their seminal theoretical work pjjj], Hurley and 
Keady showed how to get the velocity profile of the beams 
generated by such a cylinder. First, they proposed to 
consider a dimensionless parameter A = v/(2R 2 ou cot 6) 
which is proportional to the squared ratio of the oscilla- 
tory boundary layer thickness S that surrounds the cylin- 
der to its radius R. The cylinder emits the four beams of 
the well known St Andrew's configuration. Each beam 
is considered independently. Then, the longitudinal ve- 
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Ji(K) exp(-\sK 3 /R + iaK) dK, 



(2) 

where J\ is the Bessel function of the first kind of or- 
der 1 and s (resp. a) the longitudinal (resp. transversal) 
coordinate along (resp. across) the beam. Vo is the pro- 
jection of the oscillation velocity along the beam. In the 
domain S = (y /uj) 1 ^ 2 <C R, this function happens to be 
localized on the characteristics a = ±R corresponding to 
the lines of slope tan#, tangential to the cylinder. In this 
limit case, 2x4 = 8 beams are emitted from the four 
critical points of the cylinder. 

In our case, we only consider the upper characteristic 
corresponding to the single critical point C of the prob- 
lem. Hurley and Keady's expression restricted to this 
single characteristic happens to fit very well our exper- 
imental data. In the tidal experiment, A ~ 6.3 10 -4 . 
We estimate the buoyancy perturbation of the stratifi- 
cation by means of Eq. (pQ). By appropriately taking 
into account the inclination of the beam, one gets a final 
expression for the vertical density gradient that can be 
compared to our experimental data. 

Fig. [3] compares the vertical density gradient field in 
the experimental (a) and the theoretical (b) cases. The 
analytical solution of Hurley and Keady is computed on 
an imaginary cylinder tangential to the shelf break re- 
gion. The amplitude of oscillation of the cylinder corre- 
sponds to the amplitude of the tidal flow. The agreement 
is qualitatively very good. The phase opposition between 
the two parts of the beams is also well described by the 
model. 

To get a more quantitative comparison between exper- 
iment and theory, we plot on Fig. [H five vertical cross sec- 
tions (indicated in Fig. [3|) of the vertical density gradient 
field at x = —6, —3, 6 and 9 cm relative to the generation 
point. The quantitative agreement between the experi- 
mental and the theoretical profiles is very good, which 
confirms that the structure of the emitted beam only de- 
pends on the viscous boundary layer thickness S and on 
the local radius of curvature of the shelf R. 

First of all, this experiment proves that the internal 
tide generation conditions are met when the radiated 
tidal beam is aligned with the slope of the topography, 
corresponding to a critical point of emission which hap- 
pens to be amphidromic. We clearly demonstrate the 
absence of any transverse emission at that point that 
is supposed to occur in Baines analytical study Q. It 
is however important to emphasize that this disagree- 
ment might come from the absence of any sharp corner 
in the shelf break, contrary to what has been considered 
in Refs Q, @, [H[. Following this observation, we derive 
a simple analytical model derived from the theoretical 
work of Hurley and Keady [13[ and show that we can 
quantitatively estimate the density perturbation profile 
as a function of the viscosity and the local radius of cur- 



vature of the topography. 



5 F 




-10 -5 




-10 -5 5 




-10 -5 5 




-10 -5 5 

z (cm) 



FIG. 4: Vertical cross sections of the measured vertical density 
gradient (thin line) and its theoretical prediction using Hurley 
and Keady's theory (bold line). From top to bottom, the cross 
sections correspond to the abscisse x = —6, —3, 3, 6 and 9 cm. 



Since this paper focuses on the internal tide genera- 
tion, we dot not consider the evolution of the internal 
tides as they propagate further into the deep ocean. Sur- 
face reflections, in particular, are expected to generate 
internal solitary waves by creating local disturbances at 
the seasonal thermocline [H, [l5[ . 

Moreover, the emitted internal wave of tidal period 
often breaks up into internal waves of shorter period. In- 
deed, and very interestingly, a recent numerical report 
has emphasized that parametric subharmonic resonance 
may come into play in the rotating case [10[ . Experimen- 
tal work along those lines is in progress. 
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